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We extend the Lieb-Schupp theorem to the Heisenberg models with higher order interactions on 
non-frustrated or frustrated finite lattices. These lattices are constructed by even numbered rings 
with or without crossing bonds and have reflection symmetry. The results show that all ground 
states have total spin zero in wide interaction parameters region which is not covered with the 
results of the Marshall-Lieb-Mattis type arguments. 
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I. INTRODUCTION 

The Heisenberg models with higher order interactions 
have been discussed from various points of view. Hamil¬ 
tonian of the simplest model consists of bilinear and 
biquadratic exchange interaction terms. The S = 1 
bilinear-biquadratic exchange interaction model has been 
investigated in the context of the Haldane conjecture [lj 
and as models of one dimensional spin-Peierls material 
Li 2 VGe 2 0 6 J, 1 two dimensional gapless spin liquid 
material NiGa 2 S4,[l-f6] magnetism in S = 1 bosons@, 
[$] and three-flavor fermions @ trapped in optical lat¬ 
tices, magnetism in iron pnictide superconductors floUTO ] . 
and deconfined criticality and Landau forbidden phase 
transition [I3 - fl4j ■ In the case of S = 3/2 this model 
is also investigated as a model of the chromium spinel 
oxides ACr 2 04(A = Hg, Cd, Zn).[jll [l6| The S = 3/2 
bilinear-biquadratic-bicubic exchange interaction model 
is studied as a model of magnetism in S = 3/2 fermions 
trapped in optical lattices jl 7H1 and as a resource of 
the measurement based quantum computer. [20l. IfjTl In 
particular the S = 1 model was extensively studied by 
theoretical works, but models with S > 3/2 and/or bicu¬ 
bic and more higher order interactions are less studied 
and is of importance to the understanding of magnetic 
properties of the chromium spinel oxides and cold atomic 
gases in optical lattices. 

The Marshall-Lieb-Mattis theorem is one of the most 
famous exact results of quantum spin systems. In the 
case of the antiferromagnetic Heisenberg model on bipar¬ 
tite lattices with the same number of sublattice points, 
it proved that ordering energy levels, i.e., the lowest en¬ 
ergy level for allowed total spin is monotonically increas¬ 
ing function of total spin and ground state is unique 
spin singlet. (HI This theorem was extended to the case 
of the spin-iS 1 bilinear-biquadratic exchange interaction 
model. [23l - [25j] 

The Marshall-Lieb-Mattis theorem has made a last¬ 
ing contribution to check validity of a huge number of 
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results for the numerical studies of quantum spin sys¬ 
tems on bipartite lattices by using the exact diagonal- 
ization, density matrix renormalization group, quantum 
Monte Carlo simulation, etc., and now it serves as guide¬ 
lines for 1 real quantum simulators’ envisioned by Richard 
Feynman. j26|, [27J But this theorem is not applicable to 
the models on non-bipartite or frustrated lattices. In 
1999 Lieb-Schupp succeeded to prove that ground states 
of the antiferromagnetic Heisenberg model on checker¬ 
board type of the square lattice with crossing bonds have 
total spin zero. [28 h 30| Their method use reflection sym¬ 
metry of Hamiltonian, on the other hand, the Marshall- 
Lieb-Mattis theorem is based on the Perron-Frobenius 
theorem and works well if it can be find suitable uni¬ 
tary transformation which leads to same sign of off- 
diagonal matrix elements of irreducible unitarily trans¬ 
formed Hamiltonian satisfying the Perron-Frobenius the¬ 
orem. But it seems that there is no systematic method 
available to find it so far. The Lieb-Schupp theorem can 
be applied to a class of frustrated spin systems on reflec¬ 
tion symmetric lattices, but it can not give any informa¬ 
tion for the degeneracy of the ground state. 

Our purpose in the present paper is an extension of 
the Lieb-Schupp theorem to the Heisenberg models with 
higher order interactions on finite size lattices which are 
constructed by even numbered rings with or without 
crossing bonds and have reflection symmetry. As ex¬ 
plained above, the Marshall-Lieb-Mattis type argument 
does not work for non-bipartite lattices. Adding anti¬ 
ferromagnetic crossing bonds to even numbered rings in¬ 
duces a frustration of Neel order and breaks bipartiteness 
of lattices, but their reflection symmetry are preserved. 
By using this nature of lattices we will prove that all 
ground states of these models possess total spin zero in 
wide parameter region which is not covered with results 
of the Marshall-Lieb-Mattis type arguments. 

This paper is organized as follows. In section El we 
introduce some definitions and notation used through¬ 
out this paper. In section Ell to keep the paper self- 
contained, we explain a basic setup and ideas of the Lieb- 
Schupp theorem and apply this theorem to the models on 
even numbered rings to prove that all ground states have 
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total spin zero. In section El with Hamiltonian on even 
numbered rings discussed in section IIIII as a local Hamil¬ 
tonian, we construct global Hamiltonian on two dimen¬ 
sional lattices. In particular we treat square and honey¬ 
comb lattices with crossing bonds. In section 0 we sum¬ 
marize and discuss the results of sections cni and El and 
comment on the effects of the crossing bonds in infinite 
system of the S = 1 bilinear-biquadratic exchange inter¬ 
action model and physical realization of ferroquadrupole 
(spin nematic) phase in magnetic materials. 


II. DEFINITION AND NOTATION 

In this paper we study the isotropic spin-S Heisenberg 
model with up to the n-th order (1 < n <C oo) interaction 
term: 

n 

H n = - Y J2 J k(\ x -y\)( s ( x ) ■ s (y)) k , (!) 

x^y£Ak—1 

on lattice A, where — J k { \x — y |) are the fc-th order in¬ 
teraction coefficients between sites x and y £ A. The 
summation over x ^ y £ A counts every pair (once and 
once only). S(x) = (Si(x), S 2 (x), S 3 (x)) denotes spin-iS 1 
operator on site x and satisfies the usual commutation 
relations: 

[.S)(x), Sjiy'j^ — icijkSk(x')d X y- (2) 

Here we use a usual basis in which S3 (2;) is diagonalized, 
Si(x) and S3 (x) have real matrix elements and S 2 ( 2 ;) 
pure imaginary. 

This Hamiltonian can be written as the spin-S 
isotropic Hamiltonian with up to 2 n -pole interaction 
term: 

n k 

Hn = - Y Y Ik (\ x ~y\) Y °k,q( x ) 0 l,q(y)’( 3 ) 

x^y^A k= 1 q=—k 

where the Racah operators O k , q (x) (2 fc -pole operators) 
satisfy the relations: 

Oh,h{x) = ^ m )\} 1/2 ( S + (x)) k , (4) 

Oljx) = (-l)«O fc ,_,(*), (5) 

[S 3 (x),O k:q (y)] = qO k:q (x)S xy , (6) 

[S±(x),O k: g(y)\ = \Jk{k + 1) - q(q± l)O k ,g±i{x)5 xy , 

(7) 

with S±( x) = Si(x) ± iS2(x) and —k < q < — fc.[3lj| 
— /fc( \x — y |) are the 2 fe -pole interaction coefficients be¬ 
tween sites x and y. Relations between multipole inter¬ 
actions and higher powers of Heisenberg interaction are 


known to be 

1 

Y °hq( X )°i,q(.y) = s ( x ) • S (y )> ( 8 ) 

9=-l 

2 

Y °2.g(2 ’)Ol q (y) 

«=-2 

= | (S(x) ■ S(y)) 2 + ^S(x) • S(y) - ^S 2 (S + l) 2 , (9) 

and for k > 3 they are given by equations (B.20) and 
(B.21) in reference [H|. So H 2 is written as 

1 

H 2 =- Y Tl ^\ X ~y\) Y °l,Q( X )°l,q(y) 

Xy^y^A q=—l 

2 

- Y / 2(l a; -2/|) Y ° 2 ,q( X )°iq(y) ( 10 ) 

x^y£ A q—~ 2 

with 

h{\x y\) = M\x - y\ ) - (ii) 

I 2 {\ x ~y\) = ^M\x-y\), (12) 

where we have omitted a constant term in the right hand 
side of equation 110- 

Let us also introduce the total spin operator: 

^»tot ^ Qitot Qjtot ^rtot^ 

Y Si ( x )’ Y s2 ( x )> Y s ^ x n ■ ( 13 ) 

xEA x£A x£A / 

We easily see continuous symmetry of 7i n : 

[S^W.n] = [S‘ ot , H n \ = [Sf\ H n ] = 0. (14) 

III. MODELS ON EVEN NUMBERED RINGS 

In this section we discuss conditions for establishment 
of the Lieb-Schupp theorem which applies to Hamiltonian 
© on even numbered rings and properties of its ground 
state. 

Before we move forward, let us explain a setup of fi¬ 
nite size lattices A which is needed to establish the Lieb- 
Schupp theorem. Throughout this paper we consider 
A = Al U Ar which has an even number of indepen¬ 
dent sites and A can be split in two equal parts Al and 
Ar. Al and Ar are mirror images of one another about a 
symmetry plane without sites which cuts bonds between 
sites x £ Al and y £ Ar, and the collection of these sites 
is denoted by Ac if A is single even numbered rings. 

In the following we treat the models on just single even 
numbered rings, i.e., closed chains with even number of 
sites, to prepare constructions of the models on two di¬ 
mensional lattices in section Hvl 
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A. Ground state of models on even numbered rings 
and their reflection symmetry 

Following the above manner let us write Hamiltonian 
© on single even numbered rings with 2 m sites: 

= K + h% + h 2 c m , (15) 

where 

n k 

K = - E E I! x ~y\)°k, q ( x )ol iq {y), 

x^£y(z Al k= 1 q=—k 

(16) 

n k 

h R=- E E E T k{\ x - y\)° k , q (x)°l, q (y), 

k= 1 q=—k 

(17) 

n 

hl m = J2(- 1 ) k ~ 1 °k-°k\ ( 18 ) 

k =1 

with 

o L — (o L n L ••• o L 1 

a/ — i k k "> k k 15 5 _ k ) 

— ^ ( Ot-k (x) (x), @k,k—l (*t) j * * " ; Ok, — k (*t)) ; 

»GAl 

(19) 

r)R, _ r)R j 

k k 1 '“'A fc 1 j 5 _ fc ) 

— ^ ( &k(x ) {Ok,k{x ), Ok : k —1 (*t ); * i Ok, — k(x )) ■ 
i'^Ar 

( 20 ) 

^l(r) I s a collection of bonds between sites x,y £ A L (rj 
with m sites, hi™ contains parallel and crossing bonds 
between sites x £ Al and y £ Ar which have 2 m sites. 
Parallel bonds are interactions between x £ Al and 
y = x' £ Ar. Here x' means a reflection symmetric 
lattice point of x about the symmetry plane. So parallel 
bonds are perpendicular to the symmetry plane. Cross¬ 
ing bonds between a; £ Al and y ^ x' £ Ar are not per¬ 
pendicular to it. C^ (R) has 2 fc + l components which are 
given by, for -k < q < k, J 2 x e a l(r) a k{x)O k , q (x) with 
real coefficients a.k{x ) = a.k{x'). In Hamiltonian (3) in¬ 
teraction coefficients —Ik{\x — y\) depend on distance be¬ 
tween sites x and y, but, for the subsequent discussions, 
we consider Hamiltonian (1151) containing site-dependent 
interactions in hi™. 

To clarify the setup of the even numbered rings, we 
explain examples of h *™ . In the case of m = 1 /i^ ing 
is just one parallel bond, m = 2 single square with two 
parallel bonds and two crossing bonds, and to = 3 single 
hexagon with three parallel bonds and six crossing bonds. 
As one of examples, /i® ing is illustrated in FIG. [L] For 
to > 4 they are given by the same manner. 

Let us perform unitary transformation on Hamiltonian 



FIG. 1. Even numbered rings A split into equal parts Al 
and Ar which are mirror images of one another with respect 
to symmetry plane denoted by the black dashed line. As an 
example of the setup of even numbered rings, /i® illg is shown. 
It is decomposed as hi, /ir, and h%, which limitedly contains 
red, blue, and black bonds denoted by solid lines, respectively. 


(US): 

i2m _ rrt/.2m tj 
'‘'ring u “'ring 1 '' 



n k 



= K + hH - E E °IA 

(21) 


k= 1 q——k 


with 

U = exp | br E Si(%) ) , 

(22) 


V ccGAr / 


where 

we have used 



U^O\ q (x)U = (—\) k Ok,q{x), 

(23) 


for x £ Ar. The matrix elements of all matrices appear¬ 
ing in are real since Ok,±k are given by the k -th 

power of S± with real matrix elements and Ok, q are gen¬ 
erated by the repeated use of commutation relation © 
between Ok,±k and S± which have real matrix elements. 

Now we can write a ground state of /i^ g , 

= E c “/^« ® V>j8 > ( 24 ) 

a,(3 


with real coefficient matrix c a p, where {^} form a real 
orthonormal basis of S 3 eigenstates for the left subsys¬ 
tem and {^> R } are the corresponding states for the right 

subsystem. The ground state energy of h ^( g is given by 




h- 

ring 


V’) = 


TtccUi™ + Tr c'ch% 


E E Trct °t< 

k— 1 q=—k 


Rf 

9’ 


(25) 


with 


= ( 26 ) 


(°L ) a7 



E Otk{x)O k ,q{x) 
®eA L 



(27) 
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and similarly for (h]g) and . 

Following the arguments in references Fisl - fiioi . ITmL |35| 
we set c —> c\ then the right hand side of equation C5l) 
is written as 

n k 

TrctchZ 1 + W/# - ^ E Trct °M<. ( 28 ) 

k= 1 q=—k 

where we have used equation © and the cyclic property 
of trace in the third term. Here we note that Hamilto¬ 
nian m is left-right symmetric. So we can see that the 
ground state energy is unchanged and eigenstates with 
coefficient matrices c* and c + c* are also ground states. 
There exists at least one ground state with Hermite co¬ 
efficient matrix. Hermite coefficient matrix can be diag¬ 
onalized and then the third term in the right-hand side 
of equation (l?5l) are written as 

n k 

~ } ' ^ ^ ' c ll c mm |(Qfc,qr); m | i (29) 

k= 1 q=—k l,m 

in the diagonal basis of c. This expression is bounded 
below by 

n k 

-EE Emi c ™-iK°m)J 2 - (so) 

k= 1 q——k l,m 

So we can confirm that an eigenstate with positive 
semidefinite coefficient matrix |c| = \T<? is a ground state 
of h 

ui 'Sing- 


B. Singlet Ground States 

In this subsection, at first, we confirm that a ground 
state of g , i.e., ip = Uip with positive semidefinite 
coefficient matrix |c| has S'* 0 * = 0 . 

Let us introduce a tensor product of spin singlet state: 


Next, we show that all ground states of have 
S*°* = 0 even if the ground state is degenerate. Let b(x) 
be a real valued function of site x. Now we consider the 
unitarily transformed Hamiltonian under site-dependent 
field given by 


hr^b) =h? + K 


n k 

EE 

k^l q=—k 


°k, q °t q 


]ol q o* q 


9/0 



O* o- E b ( x ') 

V. x'GAr / 


1 

2 



- E b ^ 

xSAl 


( 0 &>) 


2 


1 

2 



E 

i'EAr 


b[x’) 



(34) 


Here we note that h^ g ( 0) = U^h^ g U. Following 
the argument of Kennedy-Lieb-Shastry with a trace 
inequality (3^. [36]. we get 

e ring(b) > e ring(0)> (35) 


concerning for the ground state energy of h^ g (b). It is 
required for establishment of this inequality to satisfy the 
conditions: the matrix elements of the matrix represen¬ 
tations of 0 \ q and are real and the coefficients of 
all interaction terms 0 \ q O^ are negative. 

When we choose 


b(x) 


b x £ Al, 
—b x € Ar 


(36) 


in equation ( 1 M 1 ) . it becomes 

h^{b) = h?™ - 2 mb (O z L 0 - 0, R 0 ) + 2 m 2 b 2 . (37) 

Let ip(b) be a ground state of h^ g (b). By the variational 
principle and inequality (1351) . we have 


S 

^0=0 E (-1) S_M |S, M) x ® 15, -M ) x ,, (31) 

i£A l M=—S 

where M is eigenvalues of S 3 (x) and x € Al counts every 
pair x and its reflection symmetric point x'. Ground state 
ip with |c| is written as 

^=El C la«V£®VE- (32) 

a 

We can easily see 

(ip 0 \ip) = (ipoWl^ = Tr\c\ > 0. (33) 

Since ipo has S tot = 0 and S*°* is a good quantum num¬ 
ber, ground state ip must take S tot = 0. Thus we can find 
that there exists at least one ground state with S'* 0 * = 0. 
This result makes strong in the following argument. 


(m ^'ring (b) V’(O) 

) > (m h 2 ™ g (b) m) 

= 4™ (b) > e rSg(0)j 

(38) 

which leads to 



-mb (jp{ 0 ) \0\ fi - 

Of 0 #)Uro 2 6 2 > 0 . 

(39) 


This result is independent of value of b. In order to es¬ 
tablish this inequality for arbitrary value of b it must be 

(^(0)|O z V°mE~(0))=0. (40) 


Noting that Uip( 0) = ip( 0) is ip , then we get 


ip 


E ai(x)Oifi(x) 

i£Al 


-(-!)* y ai{x’)Oi fi (x') 

ie'EAr 



(41) 
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Setting l = 1, we can see (ip 15'| ot | ip) = 0 since the above 
equation must establish arbitrary values of oq (x) for all 
x. By the rotational invariance of h^ g as is shown in 
equation (fill) , this result also holds for S\ ot and S Jj 0 *. So 
it concludes that the all ground states have S tot = 0. 

IV. CONSTRUCTIONS OF MODELS ON 

LATTICES WITH LOCAL HAMILTONIAN ON 
EVEN NUMBERED RINGS 

In the previous section we have showed that all ground 
states of on even numbered rings have S tot = 0. In 
this section we consider constructions of global Hamil¬ 
tonian on two dimensional lattices with local Hamilto¬ 
nian . Here we suppose that whole lattices are con¬ 
structed with even numbered rings, such as square lattice, 
honeycomb lattice, 1/5-depleted lattice ('CaV.iOabjddj 
1/5-depleted lattice consists of squares and octagons. 

In subsection IIV A1 we show that all ground states of 
global Hamiltonian with site-dependent interactions pos¬ 
sesses 5 tot =0. In subsection hvb] we consider mod¬ 
els on lattices without crossing bonds as a special case 
of global Hamiltonian in subsection IIV Al and determine 
conditions realizing spatially isotropic interactions. In 
subsection IIV Cl and IIV PI as examples of lattices with 
crossing bonds in this framework, we perform the same 
procedure in subsection IIVBI in the case of square and 
honeycomb lattices with crossing bonds. 


A. Generalized lattices 

In section m we have treated on even numbered 
rings and have showed that their ground states possess 
S tot _ q It is straightforward to prove that ground 
states of global Hamiltonian on generalized lattices have 
the same result. In this subsection we shortly explain it 
as follows. 

Let us write global Hamiltonian: 

n=H L +H R + H c . (42) 

He is constructed with translated copies of local Hamil¬ 
tonian y™ h,Q on an even numbered ring with 2 m sites 
in the direction parallel to the symmetry plane. [37) i.e., 

m 

= (43) 

Ac j 

where the summation for j is taken if it is necessary and 
the same applies hereinafter. The sites belongings to Ac 
are denoted by the collection of sites which are translated 
copies of sites belonging to an even numbered ring such as 
single square and single hexagon in the direction parallel 
to the symmetry plane (see FIG.2). 

Let = h™ + fig + Y™ . "Hl(r) consists of the 
collection of translated copies of h^™ g on A L ( R ) and the 


collection of bonds in h r /p Rj on 9 A L (r). The sites belong¬ 
ing to 9 A L (r) is denoted by A L (r) D Aq (the sites belong¬ 
ing to Ac = <9 Al U <9Ar). Then, global Hamiltonian on 
A L (r) is written as 

n H K)= Y ^g+ E h w (44) 

Al(r) 9A L ( R ) 

The second term in this Hamiltonian is omitted if global 
Hamiltonian is constructed with bond sharing even num¬ 
bered rings. These operations should be realized to con¬ 
struct the global Hamiltonian on the two dimensional lat¬ 
tices A = Al U Ar. We recall that Al and Ar are equal 
parts and A has reflection symmetry about the symme¬ 
try plane. Here we use ‘generalized’ lattices in the sense 
that global Hamiltonian (l42l) is constructed with con¬ 
taining site-dependent interactions in spite that H n has 
spatially isotropic interactions. 

We can easily see that global Hamiltonian (1421) also 
holds the same results in section Hill Roughly speaking, 
main differences are that, in equations (151)1 and (ESI), 
/i™/R) and local ground state ip are replaced by "Hl(r) and 

global ground state T = Y a0 C a p^<S)^ and X) Ac 
appears in the terms of parallel and crossing bonds on 
Ac- Through the same procedure in subsection IIIIBl we 
can see similar equation with respect to Ac(= (?AlU<9Ar) 
to equation m as follows. [48j 


'F 


Y. a(x)Oi, 0 (x) 


-(-!)* Y a ( x ')°ifi(x') 

x'gOAr 



o 


(45) 


and 



Y °l,o( x ) ~ (- 1 )' Y °lA x ') 

xedAi, x'COAr 



0 . 

(46) 


Let us set l = 1 and impose a periodic boundary con¬ 
dition in the direction perpendicular to the symmetry 
plane, we conclude that all ground states of global Hamil¬ 
tonian have S tot = 0 if whole lattices can be constructed 
with translated copies of Ac- Otherwise, we need to set 
other symmetry planes and impose periodic boundary 
conditions in the directions perpendicular to those sym¬ 
metry planes. 


B. Lattices without crossing bonds 

In the case of bipartite lattices, whole lattices are 
constructed with even numbered rings without crossing 
bonds, so we can construct global Hamiltonian with 
h 2 ling (m = 1), i.e., nearest neighbor pairs (x, y) only. This 
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- 4 ( 1 ) . - 4 ( 73 ) 


. - 4 ( 72 ) . - 4 ( 2 ) 

FIG. 2. (a) Square lattice with crossing bonds, (b) Checker¬ 
board type, (c) Local Hamiltonian (148II on single square with 
crossing bonds. Black solid and red dashed lines represent 
nearest neighbor interaction —7*.(1) and second neighbor in¬ 
teraction —7fe(\/2), respectively, (d) Local Hamiltonian (1551) 
on single hexagon with crossing bonds. Black solid, orange 
dashed, and blue dashed lines represent nearest neighbor in¬ 
teraction —7fe(l), second neighbor interaction — 7*, (43), and 
third neighbor interaction — 7fc(2), respectively. Black dashed 
lines in FIGs (a)-(d) represent reflection symmetric plane. 


simplest model of Hamiltonian (l42l) is written as 

v h 2 

/ ; ' ring 

MeA 

n k 

= °kA x )°k, q (y)-( 47 ) 

(x,y)£A k—1 q=—k 

Interaction coefficients in the right hand side of this equa¬ 
tion are correspond to —4(1) and we see (—l) fc 4(l) > 0 
for each k. 


C. Square lattices with crossing bonds 

In references [H, |29|, HH j39| the antiferromagnetic 
Heisenberg models with next nearest neighbor interac¬ 
tions on the square lattice and its checkerboard type were 
discussed. In this subsection, following these previous 
studies we treat the case of Hamiltonian ©■ 

Now let us consider constructions of global Hamilto¬ 
nian on the square lattice with h^ lng (m = 2). h^ ing is 
local Hamiltonian on single square with crossing bonds 
(see FIG.3-(c)). To determine the condition realizing spa¬ 


tially isotropic interactions of the models we write 

n k 2 

^ring (- 1 ) k ~ 1 ^ a k( X ) 2 °k,q( X ) 0 l, q ( X ') 

k—1 q— — k x=l 

n k 2 

-J2J2 4(i )^2°k, q { x )°k, q ( x ') 

k—1 q——k x=l 

n k 2 

= '•;,»»+EEE i-«i) 

k=1q=—kx=l 

-(-1 ) k ~ la k(x) 2 ] Ok, q ( x )°l , q ( x ')> ( 48 ) 

where — 4 ( 1 ) are denoted by coefficients of the first 
neighbor 2 fe -pole interactions. To explain the lattice 
structures and their spatial isotropy, they are illustrated 
in FIGs 2-(a),(b),(c). Noting that the factor (—l) fc ap¬ 
pears from each term in the right hand side of this equa¬ 
tion through the unitary transformation, then it can be 
seen as a combination of hf ing and h 2 ing if interaction 
coefficients of these local Hamiltonian satisfy 

-(—l) fc 4(l) + a fe (l) 2 <0, (49) 

-(-l) fc 4(l) + a fc (2) 2 <0, (50) 

for each k and we get 

(-l) fc 4(l) > 0 , (51) 

-4(1) < -4(42) < 4(1), (52) 

with second neighbor interaction coefficients — 4 (\/ 2 ) = 
(-l) fc - 1 a fc (l)a fc ( 2 ). 

When global Hamiltonian is constructed with local 
Hamiltonian (l48l) on bond sharing squares as in FIG 2- 
(a), the results in subsection IIV Al hold if —4(1)/2 < 
4(42) < 4(l)/2 and (—l) fc /fe(l) > 0. But its checker¬ 
board type is site sharing as in FIG 2-(b). So condition 
-4(1) < 4(42)/2 < 4 ( 1 ) is replaced by -4(1) < 
4(42) < 4(1). 

D. Honeycomb lattice with crossing bonds 

In this subsection we treat models on the honeycomb 
lattice with crossing bonds which is constructed with 
translated copies of local Hamiltonian on single hexagon 
(FIG 2-(d)). Let us consider global Hamiltonian with He 
given by 

/4 he = 5Z ( h °c + he) ■ (53) 

Ac 3 =2 Ac 

In FIG. 2-(d) h% is local Hamiltonian on a rectangu¬ 
lar formed by four sites 1,1', 3, 3' except dashed orange 

bonds. So inner product of 0\ and in h q is defined 

by 

0 \ = ^2 0k(x)(Ok,k{x),Ok,k-i(x),--’ ,Ofc-fc(x)) 

a;— 1,3 


( 54 ) 
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and with f3k{x') on sites 1', 3'. 

Following the previous subsection let us write the local 
Hamiltonian: 


n k 3 

^ring ~E E (- 1 ) k ^ 1 ^ a k( X ) 20 k,q{ X ) 0 l, q ( X ') 
k—1 q= — k x=l 

n k 3 

“E E ^2 T k(\ x ~ x '\) 0 k,q{ x ) 0 l,q( x ') 

k—1 q— — kx=l 
n k 

+ h c~Yl E ( _1 ) fc_1 E ^k(x) 2 0 Kq (x)0l q {x') 

k =1 q— — k £=1,3 

= h Y - m g + JlQ 
n k 

+ E E {[-A(2)-(-l) fc - 1 a fc (2) 2 ]0 M (2)0j i ,(2 / ) 

k= 1 q— — k 

+ [- 4 ( 1 ) -{-l) k ~ l a k {xf 

£=1,3 

— (—l) fe_1 /3fc(^) 2 ] °k,q{ x ) 0 k, q ( x ')} . (55) 

where —4(1), — 4(v / 3), and —4(2) are denoted by first, 
second, and third neighbor interaction coefficients of 2 fc - 
pole interactions, respectively, as illustrated in FIG. 2- 
(d). Similar to the previous subsection, right hand side 
of this equation can be seen as a combination of h® ing and 
Hq along with h 2 ing and its analogue if 

— (—l) fe 4(2) + ttfc(2) 2 < 0, (56) 

-(-l) fe 4(l) + a k ( l) 2 + p k (l ) 2 < 0, (57) 

-(-l) fe 4(l) + a fe (3) 2 + /3fe(3) 2 < 0, (58) 

and conditions on interaction parameters satisfying spa¬ 
tial isotropy in A are given by 

-4(v / 3) = (-l) fc - 1 Q fe (l)a fe (2) = (-l) fc - 1 a fc (2)a fc (3), 

(59) 

-4(2) = (—l) fe_1 [a fc (l)a*(3) + &(1)&(3)], (60) 


«fc( 1) = a k ( 3), and /4(1) = Pk( 3) for each k. From these 
equations and inequalities we get 


1 ^ 4(2) ^ i l ( h(Vs) \ 

2 - 4(1) - 4 + \| 16 + ^ 4(1) ) ’ 

(_l) fc 4£) > ( _ 1)fe4(2) > o, 


(61) 

(62) 


where we have replaced 4 ( 1 ) by 4 (l )/2 since global 
Hamiltonian is constructed with local Hamiltonian (15511 
on bond sharing hexagons. 


V. SUMMARY AND DISCUSSIONS 

We have discussed the Heisenberg models with higher 
order interactions or multipole interactions on finite lat¬ 
tices with reflection symmetry written as in the form of 


Hamiltonian (usd or m and have found that there ex¬ 
ists at least one ground state with 5' tot = 0. Moreover 
imposing a periodic boundary condition in the direction 
perpendicular to the symmetry plane, we have confirmed 
that the all ground states possess S tot = 0 even if the 
ground state is degenerate. These results are a straight¬ 
forward extension of the Lieb-Schupp theorem to these 
models. 

For establishment of the results in subsections IIII Bl 
and IIV Al we did not put any restrictions on signs or val¬ 
ues of interaction coefficients — I k {\x — y\) of /i™ and h 
except their reflection symmetry (see Hamiltonian 03 - 
CiIl) and (l42l) - (l44l) 'l. These coefficients are not essential 
to our results. In subsections IIV Bl IIV Cl and IIV Dl they 
were determined by the conditions in order that pos¬ 
sesses spatially isotropic interactions as in FIG. 2. On 
the other hand restrictions on interaction coefficients of 
parallel bonds and crossing bonds in hQ n come from the 
establishment of inequality (135(1 and similar inequalities 
for the models in this paper. |48( Therefore we have no 
idea for improvement of these restrictions. 

In this section, we summarize and discuss the results 
in section mu and m which are divided into the models 
on lattices with and without crossing bonds. 


A. Lattices without crossing bonds 

In this case, models are constructed with local Hamil¬ 
tonian h 2 ing . Typical examples of the whole lattice are 
bipartite lattices such as hypercubic lattice, honeycomb 
lattice, and 1/5-depleted lattice. In the case of the 
1/5-depleted lattice we should set the symmetry plane 
which intersects octagons and squares. The results hold 
if (—l) fe 4 (l) > 0 for each k. 

In the following we shall explain comparisons with 
the results of the Marshall-Lieb-Mattis type argument. 
Hi = Hi with — 4(1) = —4(1) > 0 is the spin-5 1 antifer¬ 
romagnetic Heisenberg model. The Marshall-Lieb-Mattis 
theorem assures that its ground state is unique and has 
5 tot = 0. So the result given by the Lieb-Schupp the¬ 
orem is completely covered by the Marshall-Lieb-Mattis 
theorem with uniqueness of the ground state. 

74 is equivalent to bilinear — 4 ( 1 )(= — 4 ( 1 ) — 

34(l)/4) and biquadratic — 4(1)(= — 3 l 2 (l)/ 2 ) ex¬ 
change interaction model H 2 ■ From the previous studies 
of the Marshall-Lieb-Mattis type argument it was known 
that the same results hold for 74 with 5 = 1 in the region 
4(1) > 4(1), 4(1) > 0(4(1) 7 ^ 0 ) and with 5 > 1 in 
the region 0 < 4(1) < —4(l)/25(5— l). [23l - f25j On the 
other hand, our results based on the Lieb-Schupp the¬ 
orem show that all ground states have 5 tot = 0 in the 
region 4(1) > 24(1), 4(1) > 0 for any 5. For 5=1 
purely biquadratic interaction model ( 4 ( 1 ) = 0 , 4 ( 1 ) > 
0) satisfies SU(3) symmetry and its ground state is de¬ 
generate. Our results can conclude that all degenerate 
ground states possess 5 tot = 0. In the case of 5 > 1 our 
results extend the region which one can conclude ground 















states with S'* 0 * = 0. These results are summarized in 
FIG. 3. 

For 1-L n with n > 2, our results also hold if 

(-l) fe 4(l) > 0 for each k. As far as we know, the re¬ 
sults of the Marsliall-Lieb-Mattis type argument does not 
exist. 

Our study is concerned with models on finite lattices 
with reflection symmetry and their ground states pos¬ 
sess S*°* = 0, but in infinite volume limit continuous 
symmetry breaking may occur. The antiferromagnetic 
Heisenberg model on bipartite lattices in two or more di¬ 
mensions is known to exhibit Neel long range order in 
its ground state. For the d-dimensional hypercubic lat¬ 
tice, it was rigorously proved in d > 3 for any 5 and 
in d = 2 for S > lf34, iH EuJ, and for the honeycomb 
lattice for S > 3/2 [lj - Ground state phase diagram of 
5 = 1 bilinear-biquadratic model (if 2 ) with J 2 (l) > 0 
on the square or the simple cubic lattice is considered as 
follows- SM The region J 2 (l) > Ji(l) > 0 is the fer- 
roquadrupole (spin nematic) phase, Ji(l) < 0 the Neel 
ordered phase, and 0 < J 2 (l) < J 2 (l) the ferromagnetic 
phase, which are illustrated in FIG. 3-(a). In our param¬ 
eter region there exist Neel long range order and ferro- 
quadrupole long ranger order in infinite volume ground 
state, which were rigorously proved in parts of these pa¬ 
rameter region. JH, [4f|[47| These proofs were given by the 
method of infrared bounds whose key inequality is upper 
bound on the Fourier transformed correlation function in 
whole momentum space which is derived from inequality 
(13511 and similar ones for the ground state energy or anal¬ 
ogous inequality for the partition function. 




FIG. 3. (a) S = 1 and (b) S > 1 ground state phase di¬ 
agram of the bilinear-biquadratic exchange Hamiltonian on 
finite lattices without crossing bonds. Our results based on 
the Lieb-Schupp theorem mean that all ground states have 
5 tot = 0 in the dark gray region. The results of the Marshall- 
Lieb-Mattis type arguments mean that unique ground state 
has 5 tot = 0 in the light gray region. Three colored regions 
in inner semicircle in FIG. (a) are expected phase diagram 
of d, > 2 dimensional hyper cubic lattice in infinite volume 
limit. Pink, orange, and blue region is ferromagnetic (FM), 
ferroquadrupole (FQ), and Neel ordered (AFM) phase, re¬ 
spectively. 


B. Lattices with crossing bonds 

In this case, models are constructed with translated 
copies of h^ g with m > 1 on rings. These rings have 
crossing bonds. In subsection IIV Al we have proved that 
global Hamiltonian (l42l) holds the same results in sub¬ 
section ME These models have site-dependent inter¬ 
actions. So in subsections IIV Cl and IIV Dl we consid¬ 
ered conditions which Hamiltonian (H^l) possesses spa¬ 
tially isotropic interactions on the square and the hon¬ 
eycomb lattices with crossing bonds as in FIG. 2. For 
the non-checkerboard square lattice it is realized in the 
region —/*,(l )/2 < Ik(V 2 ) < /&(l )/2 and checkerboard 
type —Jfc(l) < Ik(V 2) < Jfc( 1) for each k. The result for 
the honeycomb lattice is given by inequalities m is, 
and setting |/fc(l)| = 1 it is illustrated in FIG. 4. 

Only about the frustrated antiferromagnetic Heisen¬ 
berg model on the square and honeycomb lattice with 
crossing bonds, we explain relation between our results 
on finite lattices and the results of various theoreti¬ 
cal studies on infinite lattices. There exist detailed 
reviews of these models on square lattices in the pa¬ 
per [49J (see also references therein). In infinite vol¬ 
ume limit, within our parameter region, the ground 
state phase diagram expected to be valid is summa¬ 


rized as follows. For the non-checkerboard (checker¬ 
board) square lattice, in the case of 5 = 1/2, the region 
0(0) < Ji(\/2)/Ji(l) < 0.4(0. 8 ) is the Neel ordered phase 
and 0.4(0. 8 ) < Ji(y/2)/ Ji(l) < 0.5(1) the quantum para¬ 
magnetic phase without magnetic long range order, and 
in 5 = 1, 0(0) < Ji(y/2)/Ji(l) < 0.5(1) the Neel or¬ 
dered phase. In the case of 5 = 1 and non-checkerboard 
type, by using the method of infrared bounds, the ex¬ 
istence of Neel long range order was rigorously proved 
in small J 1 (V2)/Ji flV [38(| Ground state phase diagram 
of 5 = 1/2 on the honeycomb lattice (non-checkerboard 
type) was obtained in reference 0, [5l| . Our parame¬ 
ter region also seems to be contained in the Neel ordered 
phase and the quantum paramagnetic phase. 

In this paper we focus on one dimensional rings and 
two dimensional lattices. Application of these results to 
three dimensional lattices can be easily extended. In that 
case, as a simplest example, we can consider local Hamil¬ 
tonian on single cubes with crossing bonds and it should 
be written as ft® ox . 

Lieb-Schupp called FIG. 2-(b) pyrochlore checkerboard 
since it is a two dimensional projection of a three dimen¬ 
sional pyrochlore lattice. But this framework is not ap- 
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\h{2)\ 



FIG. 4. Ground states of T-L n on honeycomb lattice with 
crossing bonds have S tot = 0 in the gray parameter region 
for each k which is determined by inequalities © and EH- 
Jfc(l) is scaled as —1 if k is odd, otherwise 7^(1) = 1, and 
— 1/2 < 7&(2) < 0 if k is odd, otherwise 0 < Ik{ 2) < 1/2. 


plicable to the pyrochlore lattice unfortunately, because 
it lacks reflection symmetry. They also called equation 
(14111 with l = 1 quantum analogue of ice rule in the con¬ 
text of the correspondence between Ising like ferromag- 
net with crystal field anisotropy on the pyrochlore lattice 
and configuration of four hydrogen atoms around an oxy¬ 
gen atom in ice. (28l . [29t | In that point of view equations 
HUD and (1461) are generalization of ice rule to any even 
numbered frustrated units and 2 2i_1 pole moment higher 
than dipole. 

In the following we shall comment on effects of crossing 
bonds on stability of the Neel ordered phase and the fer- 
roquadrupole (spin nematic) phase. As was illustrated 
in FIG. 3-(a), there exist the Neel ordered phase and 
the ferroquadrupole phase which are separated by the 
line at Ji(l) = 0. By adding antiferromagnetic crossing 
bonds to the square lattice, it is clear that the Neel or¬ 
der exhibiting anti-alignment of spin is not stable. On 
the other hand, the quadrupole order is not alignment 
of spin but nematicity of spin, and from equation (1231) 
it can be seen that Ok,o with even k is even parity with 
respect to time reversal. The ferroquadrupole order is 
uniformly aligned nematic and does not seem to be suf¬ 
fer from geometrical frustration. So stability of the fer¬ 
roquadrupole phase is not affected by frustration due 
to antiferromagnetic crossing bonds unlike the Neel or¬ 
dered phase. Now we set next nearest neighbor inter¬ 


actions — Ji(\/2) = — and —J 2 (\/2) = —aJ 2 ( 1) 

with 0 < a < 1/2, then FIG. 3-(a) is expected to 
be changed as follows. Phase boundary J 2 (l) = Ji(l) 
is unchanged by adding crossing bonds with ferromag¬ 
netic and ferroquadrupole interactions since the region 
0 < =72(1) < Ji(l) is saturated ferromagnetic ground 
state. On the other hand, phase boundary Ji(l) = 0 
closes to the antiferromagnetic Heisenberg point as a ap¬ 
proaches 1/2, i.e., by adding crossing bonds with anti¬ 
ferromagnetic and ferroquadrupole interactions, the fer¬ 
roquadrupole phase becomes dominant and the Neel or¬ 
dered phase is suppressed if there do not exist differ¬ 
ent phases between these two phases. As for stability 
of the Neel order, the triangular and pyrochlore lattices 
are slightly different situation from the square lattice 
with antiferromagnetic crossing bonds. The ground state 
phase diagrams of S' = 1 bilinear-biquadratic exchange 
model on the triangular and pyrochlore lattices are ob¬ 
tained in reference [U, mm and the same situation in 
the above scenario is shown. 

Finally we shall propose the physical realization of the 
ferroquadrupole phase in magnetic materials. Usually 
biquadratic interaction is small as compared with bi¬ 
linear interaction and the ferroquadrupole phase is un¬ 
physical in magnetic materials. In reference [3] Mila and 
Zhang proposed a mechanism leading to a significant bi¬ 
quadratic interaction in S = 1 systems as follows. The 
super exchange interaction between atoms with three 
orbitals and two outer electrons per atom, which con¬ 
sists of the two singly occupied doubly degenerate or¬ 
bitals with the lowest energy and an unoccupied orbital 
with slightly higher energy. The virtual electron tran¬ 
sition via the higher energy orbital favors ferromagnetic 
spin interaction, which compensates largely the antifer¬ 
romagnetic superexchange interaction. As a result, the 
biquadratic interaction becomes predominant relatively. 
Thus we expect that highly frustrated antiferromagnetic 
materials with biquadratic exchange interactions origi¬ 
nated from the Mila-Zhang mechanism may exhibit the 
ferroquadrupole phase. 
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